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Abstract 

The quantum fluctuations of the quark condensate are studied in a 
Nambu Jona-Lasinio model. Two Lorenz invariant regularizations are 
considered: a sharp 4-momentum cut-off and a soft gaussian regulator. 
The quantum fluctuations of the quark condensate are found to be 
large although chiral symmetry is not restored. Instabilities of the 
ground state appear when the system is probed by a source term 
proportional to the squared quark condensate. The instabilities are 
traced to unphysical poles introduced by the regulator and their effect 
is greatly enhanced when a sharp cut-off is used. 



1 Introduction. 



We study the effect of quantum fluctuations of the quark condensate on the 
physical vacuum. We use an SU (2)^ Nambu Jona-Lasinio model in the chiral 
limit. Two ultra-violet regularizations are considered: 

• regularization using a sharp cut-off, in which the quark propagators are 
set to zero when > A 2 , where k^ is the euclidean 4-momentum of 
the quark propagator; 
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• regularization using a smooth gaussian regulator which will be de- 
scribed below. 

We introduce a current quark mass m as a source term mifjip to calculate 
the quark condensate ($ip^- We introduce another source term ~j (^T a i/^j 

which allows us to calculate the expectation value ^ (^T a ij?j \ of the squared 

quark condensate. We shall however maintain j at a finite value in order to 
study the response of the system to a deformation induced by this source 
term. We include the 1/N C effects due to the quantum fluctuations of the 
meson fields. In the quark language, this means that we include the exchange 
(Fock) term as well as the (RPA) ring diagrams. 

Although the quantum fluctuations of the meson fields do not restore 
chiral symmetry, we do find surprisingly large fluctuations of the quark con- 
densate. We also find that the effective potential is very sensitive to the 
shape of the regulator. Apparent instabilities appear which we show to be 
artefacts sharp cut-of used in conjunction with the relatively low values of 
the cut-offs, used in chiral quark models. 

The quantum fluctuations of the meson fields in the Nambu Jona-Lasinio 
model have been studied for several years |J, @, @, [0- In these studies, 
the quark and meson loops were regularized with different cut-offs. Since 
both the quark and meson loops diverge, their relative contribution could be 
adjusted at will by a proper choice of the cut-offs thereby making it impossible 
to estimate the importance of the quantum fluctuations of the meson fields. 

In this study we consider the Nambu Jona-Lasinio model to be a quark 
model in the sense that all physical processes can be expressed in term of 
Feynman graphs involving only quark propagators. The meson fields which 
are introduced in the process of bosonization are mere intermediate quantities 
introduced to calculate the partition function. When the quark propagators 
are regularized, a single cut-off regularizes both the quark and meson loops. 
No further regularization is required for higher order loops. This approach 
has also been adopted in Ref. for example, with the exception that a 
3-momentum cut-off was used. We shall show that results obtained with 
3 and 4-momentum cut-offs can differ even qualitatively. We show that 
regularization with a 4-momentum cut-off introduces a non-locality which 
makes the ground state energy unbounded from below, or, equivalently, which 
makes the model acausal. We relate this instability to the unphysical poles 
of the quark propagator which are introduced by the regulator. 
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We consider the regularization of a model to be a physical phenomenon 
and not simply a way to be rid of infinities wherever they turn up. Because 
of this we regularize the model action from the outset before calculating the 
loop integrals and we shall see that this makes a significant difference with 
the common practice of regularizing the infinities of the loop integrals which 
are deduced from an unregularized action. We derive all quantities in terms 
of a regulator which is a function of the squared euclidean 4- momentum k 2 . 

The paper is composed of three parts. In the first part we set the notation, 
we explain how the calculation was performed and we define the relevant 
range of the model parameters. The second part is devoted to a discussion 
of the quantum fluctuations of the quark condensate. The last part discusses 
the instabilities which are displayed by the response of the system, subject 
to constraint proportional to the squared quark condensate. 



2 The model euclidean action. 

The model is defined by the euclidean action: 

I m (q,q) = {q~\-idnln + rmr\q) - — J d 4 x ((q \r\ x) T a (x \r\ q)f (1) 

It involves a quark field q(x) = (x\q). The euclidean Dirac matrices are 
7^ = 7^ = (i/3,7). The model includes a regulator r which is assumed to 
be diagonal in /c-space: (k \r\ k') = Skk'r(k). We consider both a gaussian 
regulator: 

k 2 

r(k) = e~^? (2) 

and a sharp cut-off: 

r (Jfe) = 1 if k 2 < A 2 r (k) =0 if k 2 > A 2 (3) 

The use of the the sharp cut-off is tantamount to the calculation of Feynman 
graphs in which the quark propagators are set to zero when their euclidean 
4- momentum k 2 exceeds the cut-off A 2 . The bra-ket notation for the Dirac 
fields is: 

(q + m\q)= J d 4 xq (x) (-id^ + m) q (x) 

(x\r\q) = d 4 y (x\r\y) q(y) (q \r\ x) = d 4 y q (y) (y |r| x) (4) 
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We use the current quark mass m to evaluate a regularized quark condensate 
and that is why the current quark mass m is multiplied by the regulator. The 
T a = (1, i^f^f) are the generators of chiral rotations. The number of flavors is 
Nf and there are Nj generators T a . We assumed that the coupling constant 

2 

jj- is inversely proportional to N c . 

The partition function W is given by the euclidean path integral^ 

-W(m) = f n(n\ n (n\ 



D(q)D(q)e- 1(q > q) (5) 

At zero temperature and for an infinite translationally invariant system, 
W = Qe where Q is the space-time volume Q — J d^x 1 and where e is the 
energy density e = y, which is the energy per unit volume of the system in 
its ground state. 

The model is regularized in the same way as the effective quark model 
which has been derived from a study of the propagation of quarks in an 
instanton liquid * n which case both the shape of the regulator and the 

value of the cut-off are derived. However, the model action (|1|) we are using 
is not exactly the same as that derived from the instanton liquid ||. The 
model described by the action ([j]) has been actively investigated in both the 



soliton M and meson sectors ]T0] O, [12 



We now adopt a notation which adds considerable transparency to the 
manipulations made below. We define an interaction V by its matrix element: 

(xa \V\ yb) = (yb\V\ xq) = -S ab S (x - y) ^— (xa V' 1 yb) = -5 ab 5 (x-y)^ 

(6) 

We define a delocalized quark field tp (x) and the corresponding bilinear forms 

ip(x) = J d^x (x\r\y) q(y) $ (x) T a ip (x) = J d^y d^z q (y) (y |r| x)T a (x \r\ z) q (z) 

(7) 

In this notation, the action (fl|) takes the form: 

im (q, q) = (ql-id^in + rmr \ q) + \ (^r^) v {^pv^ (8) 

The chiral limit is defined to be m — ► 0. 



1 Pedantically, the partition function is Z = e w . In the following we shall call W the 
partition function without fear of confusion. 
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3 The constrained system and the effective 
potential. 

The current quark mass m serves as a source term to calculate the normalized 
quark condensate ($ip\ In order to calculate the quantum fluctuations of 
the quark condensate, we introduce an extra source term: 



(Vr a ipy = -j j d A x (ipipy + (^i"/sT a ip) 



(9) 



which acts as a constraint on the system and which allows us to calculate 
the expectation value (^(^T a ifj \ of the squared condensate. We prefer to 

work with the squared condensate (tJjT a tp^J rather than with (ipip) because 
ipip is only one of the four components of the chiral 4- vector ipT a %p. If the 
system chooses to vibrate in the direction ipip defined by the ground state, 
it will do so because the amplitudes of the vibrations in the four directions 
are independent. But there is no reason to prevent the system to vibrate in 
the other three directions. 

The constraint is introduced into the action (§) and we define the con- 
strained system by the partition function: 

e -W(j,m) = J D ( q ) D (g) e -/m(g,<2)+§(^)i(^rv>) = J D ^ D (^) e --0.mfo,3) 

(10) 

where: 

J i,m {q, q) = (g |-«Vfa + rmr \ <?) + \ (V' r V') {V - j) (11) 

is the action of the constrained system. 

The ground state expectation value / (^jjT a ipj ^> of the squared condensate 
in the constrained system is given by: 

1 //- \2\ dW (j,m) 



-^((V5r v 



2 \V ^) I dj 
and the quark condensate is given by the expression: 

dW (j, m 



(12) 

m=0 



(13) 

m=0 
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where Q is the space-time volume in euclidean space. 

The constrained system, described by the partition function W (j, m) is 
not the same as the system described by the partition function W (m) because 

This is why, in the presence of the constraint, the energy density of the 
system is defined in terms of the effective potential: 



W (j, m) + ~ j ( (#^) 2 ) = W (j, m) - j 8W U m) 



dj 



(14) 



m=0 



Viewed as a function of j (for a fixed m), a stationary point of the action 
occurs when: 

dW (j, m) . d 2 W(j,m) 

dj ~ ~ 3 W 2 ( j 

so that the effective action is stationary when j = 0, that is, in the absence 
of the constraint. This is true whatever approximation we use to calculate 
W(j,m). 

The effective potential allows us to map out the energy of the system 

as it deforms under the effect of the constraint \j NjT a il^j . The choice of 
the constraint used to probe the energy surface of the system is, of course, 
arbitrary. Our choice is justified by the fact that, as we shall see in section 
pj| the system offers a relatively soft response to the constraint and in some 
cases it actually displays an instability of the ground state. 

4 Bosonization in terms of local fields. 

Bosonization is simply a convenient way to calculate the partition function 
W (j, m). We introduce local auxiliary fields (p a (x) and we consider the new 
euclidean action: 

Ij,rn{q,q,¥) = (q\~id^^ + rmr + np a T a r\q) - -if (V - j)' 1 <p (16) 
Consider the partition function W' (j, m) defined by the path integral: 

e -W'(j,m) = J D D ^ D e -Ii, m (q,q, V ) ( 17 ) 
The integral over ip yields: 
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where tr denotes a trace in (x, a) space: 

trl = nN 2 f (19) 

(because there are Nj generators r a ). From ([l?]) and (|Tj|) we deduce a 
relation between the partition functions W (j, m) and W (j, m) : 

W (j, m) = W (j, m) + -tr In (V - j) (20) 

If we integrate out the quarks in the path integral (|T7|), we get the partition 
function W (j, m) in the form: 

e -W(j,m) = J D ^ e -IjMv)-^trHV-3) (21) 

where Ij >m (tp) is the so-called "bosonized action": 

Ij,m (f) = ~Tr In [-id^ + rmr + np a Y a r) - ~tp(V - if (22) 

The trace Tr is over the variables (space-time x, Dirac indices, flavor and 
color) which define the quark field: 

Tr 1 = 4N c N f n (23) 

It is convenient to shift the constituent quark mass to the interaction term 
by writing m + ip a T a = (p' a T a - The bosonized action ( p2|) becomes: 

Ij,m (<f) = -Tr In i-id^ + r Va T a r) -~(<p-m)(V- j f 1 (<p - m) (24) 

where dropped the prime on (p. The partition function is then given by the 
expression fl2~i~D in which Ij >m ((p) is the action (|24]) . In the expression (^4|), 
m stands for the vector m a with components m a = (m, 0, 0, 0). 

It is usual to calculate the effective potential from the texbook formula 
Q: 

^^'-('^"W (25) 

obtained from a Legendre transform which relates a source term to the ex- 
pectation value of the field. We do not use this formalism because, in our 
case, the operator & q™q^ has negative eigenvalues. This is easily seen by 
considering cuts in the Mexican hat shaped I m (if) . 
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5 The classical approximation to the constrained 
system. 

The classical approximation consists in approximating the path integral (|2T|) 

as: 

W(j,Tn)=Ij,mM (26) 

where <pj is a stationary point of the action Jj >m (<p), defined in (p^). In the 
classical approximation, we also neglect the constant |tr In (V — j) which we 
will include in section ^| together with the contribution of the field fluctua- 
tions. The term "classical" will be used throughout although, at the quark 
level, it includes the quark loop. As it turns out, it is the classical approx- 
imation to the constrained system determines the response of the system to 
the constraint and the quantum fluctuations of the fields are only small cor- 
rections. The classical approximation is the leading order contribution in N c 
and, in the quark representation, it corresponds to the Hartree approxima- 
tion. 



5.1 The gap equation and the relation between j and 
M. 

A stationary point of the action (^|) occurs for ip a = (M, 0, 0, 0) where M is 
the solution of the so-called gap equation: 

Mir 4 M 
M — m 2\l —of l +r qL M' L M—m 

where qm = 9m (q = 0) and #m (q) is the function QB.6|) , defined in appendix 
p. We shall denote by Ij im (M) the action at the stationary point. An explicit 
expression is given in appendix [A]. 

Let Mo be the solution of the gap equation in the absence of a constraint 
(j = 0): 

= - 4AyV 'lCT"« < 28) 

This equation relates M to the interaction strength V. The minus sign 
indicates that the interaction V is attractive. Otherwise chiral symmetry 
would not be spontaneously broken and ip a = would be the only stationary 
point of the action. 
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We can calculate j from the equation: 

j = L_) (29) 

4N c N f \g M 9mJ 4:N c Nj \Mg M M g M J 

In practice, we start by choosing a value of M , which determines V. 
We then choose a value of M which determines j. We can easily derive the 
relation: 

which shows that M is a monotonically increasing function of j so that it 
makes little difference if we plot the effective potential as a function of j or 
M. The choice of M is more transparent. 

5.2 The classical estimates of (i/ji/j) and ^(^r a^^f^j in 
the chiral limit. 



In the classical approximation, the quark condensate (|13"D is given by the 
expression: 



- 9/ "» (M) 



= Vl(V - M = —4QN c NfMg M (31) 

m=0 



cZass (9?7l 

The classical approximation to the squared quark condensate (|T2|) is: 



niliprM ) = — ^ — n M 2 {w c N f g M y (32) 

where we used the fact that, for a given value of j, M is a stationary point 
of the action Ij m (M) . 

It follows that, in the classical approximation, the fluctuation of the quark 
condensate is zero: 

(Uv^f) =(W) 2 (33) 

\ V J I class X 1 class 

as expected. The quantum fluctuations of the quark condensate are due to 
the quantum fluctuations of the fields and they are introduced in section 
§ 

The gap equation fl27|) relates M to the classical quark condensate and 
M is simply a Hartree insertion: 
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M =(V-j) (^V) 




(34) 



class 



5.3 The classical effective potential in the chiral limit. 

Because M is a stationary point of the action Ij m (M), the classical approx- 
imation to the effective potential (|T4|) is: 



T d ass = I j>m (M) - /W M ) = j. m ( M ) + jifi (V - J)' 2 M 2 (35) 



Using (|27| ) and (p9|) , the classical effective potential can also be expressed as 



the following function of M: 

T d ass (M) = I hm (M) + ifi (4N c Nf) M 2 g M (l - (m = 0) (36) 

2 \ 9 M ) 

An explicit expression for (M) is given in appendix |A[ 



6 Inclusion of the quantum fluctuations of 
the fields if. 

A saddle point evaluation of the partition function yields: 
W(j,m) = I hm (< P3 )+ 1 -trln 6%m{ip) 



+ ~tr\n(V-j) (37) 

<p=<Pj 



2 5ip5ip 

where tpj is the stationary point of Ij jTn (ip), which is determined by the gap 



equation (|27|). 



6.1 Calculation of W (j, m). 

To calculate the second term of ([37]), we need to evaluate the inverse meson 
propagator matrix: 

K *&v) = T^rfrn i ( 38 ) 

dLf a (x) btfb (y) 
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at the stationary point of (if). From the action (H), we see that: 

K-^U-iV-j)- 1 (39) 
where II is the polarization function (Lindhardt function): 

S 2 

n a6 {x, y) = -- — — — — -r-zTr In {-id^ + rmr + np a Y a r) 
5(p a (x) 6(p b (y) 

= Tr — — ■ — - — ■ — \x) T a (x\ -— ■ — - — ■ — \y) T b (y\ 

-io^Ih + rmr + np a T a r -«<V»t + rmr + rip a T a r 

(40) 

The partition function fl37|) becomes: 

(41) 



W (j, m) = I j>m fa) + l -tr In (1 - II (V - j)) 



The matrices II and V, and therefore K, are diagonal in momentum space 
and in the flavor indices: 

qa K~ l q'bj = 5 a bS qq 'K^ 1 (q) (qa \U\ q'b) = 8 ab S qq 'U a (q) 

In particular: 

n a=0 (?) = U S (q) = AN c Nj ^-q 2 f 2 (q) + M 2 fj» (q) + f£ (q)) - g M (q) 

n a=1)2 , 3 (?) = Il P (q) = 4N c Nf (±q 2 f 2 R 2 (q) + M 2 (/» (q) - f£ (?)) - g M (q) 

(42) 

where the functions f^f (q) and gM (?) are defined by the expressions ( |B.5| ) 
and (p.6| ) of appendix |B[ The regulator r ensures that the polarization 
functions II (q) vanish when q» 2A. 

If we use the gap equation ( p7]) to express (V — j)" 1 in terms of gM = 
9m (? — 0), we obtain analogous expressions for the inverse meson propaga- 
tors: 

K~} (q) = K s ' (q) = 4N c Nj Qg 2 ffi (q) + M 2 (/* (q) + f£ (q)) - g M (?) + J^9m 

Kll,2,s (?) = Kp 1 (?) = WJf, (^q 2 fH (?) + M 2 {f» (q) - f$ (q)) - g M (q) + jJ^9m 

(43) 
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The pion remains a Goldstone boson even in the constrained system be- 
cause, in the chiral limit, K~l x 23 (q) — > 0. This is an important feature of 

the choice we have made of the constraint which does not break the chiral 
symmetry of the lagrangian. 

The partition function (O) is then: 



W(j,m) = I jjm (M 2 ) 
+~ X> (l + fafM (?) + M 2 (ft (q) + flj (?)) ~ 9M (q) 



Nj-1 
4 " ; 



_ (44) 

It is because we took the trouble to keep track of the term \tr In (V — j) that 
the sums over q converge. Another way to write W (j, m) is: 

W (j, m) = I 3 m (M) + 1 V In (45) 

6.2 The Fock term (exchange energy) and the ring di- 
agrams. 

Consider the expression ([IT]) of W (j, m) . Using the explicit form (^) of 
the action Ij tm (M) as well as the relation (0), we can write the partition 
function W (j, m), in the chiral limit, as follows: 

W (j, m) = -Tr In (-id^ + r 2 M)~ (^) {V - j) ($i/>)+ -tr In (1 — II (V — j)) 

(46) 

We can express the second term, together with the expansion of the loga- 
rithm, in terms of Feynman graphs. The second term of ( ^6[ ) is the direct 
(Hartree) energy: 



<p=ipj 



{V-j)($i/>) 





(47) 



The expansion of the logarithm in ( |4*6| ) generates the exchange (Fock) energy 
as well as the ring diagrams: 
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-trln(l-U(V-j)) 



~trU (V - j)+]trU (V - j) U(V - j)~trU (V - j) U (V - j) U (V - 
2 4 o 




+ 




+ 



(48) 

The terms of (^) are the next to leading order in N c . The first term, which is 
the exchange (Fock) term, is special in that it does not induce any correlations 
in the quark wavefunction. It is easy to see that the exchange term is more 
sensitive than the ring diagrams to the high momenta running in the meson 
loop. As a result, when a sharp cut-off is used, it is the exchange term and 
not the ring diagrams which dominates the next to leading order contribution 



+ 



6.3 Expressions for the quark condensates and the quark 
condensates in the chiral limit. 



The ground state expectation value of the squared condensate (|I~2| ) can be 
calculated from the expression: 



l n // 7r a 2 \ _ d km(M) 



dj 



m .„ Sj SAP 4N e N,g M ) 1 ' 



where is given by (p0|). The contribution of the exchange term to the 
squared condensate is: 

n^riyfM>-^ nM ) <5o) 
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The quark condensate flT3"|) can be calculated with the help of (f4"T|): 

(51) 



n (^) = = A (j,, m (M) + ifr in (i - n (v - j)) 



M 



We calculate the derivative avt ^ m - > keeping j constant. However, we must 
remember that M depends on both j and m and therefore a contribution 
arises from the change in M when m is varied because, in contrast to Ij. m , 
W (j, m) is not stationary with respect to variations of M . When m — > 
m + 5m, (V — j) remains constant and: 

5-tr \n(l - (V-j)U) = --tr V ~ 3 .. TJ 5U= -trKSTl = -^trK- 6U 



2 v v ' 2 1 — (V — j) II 2 2 5m 5M 2 

(52) 

The quark condensate is thus: 

From the gap equation Q27|) we see that ^ = (jjffs)" 1 - 

The contribution of the exchange term to the quark condensate is: 

\^^/exch y J) 5m 2tt 5M 2 1 ; 

The effective potential (|TJ]) can be calculated from the expression 

r^(j>)+i^(«) 2 ) (55) 

where j is given by (|29|) and | <^ [iJ)Y a il?j \ by ( |49|) . 

7 The relevant range of parameters. 

We consider now the parameters of the model. We calculate all quantities 
in units of the cut-off A which appears in the regulators (|2]) and (|3|). This 
is convenient because, for example, the effective potential is proportional to 
A 4 but otherwise it depends only on the ratios M /A and M/A. Thus, the 
stability of the system depends on the single parameter M /A and, of course 
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Mo 
A 


Z s 


Z P 


U (MeV) 
forM = 300 MeV 


U (MeV) 
forM = 400 MeV 


0.2 


0.171 


0.189 


130 


173 


0.4 


0.0825 


0.0964 


93.2 


124 


0.6 


0.0433 


0.0533 


69.7 


92.4 


0.8 


0.0239 


0.0307 


52.6 


70.1 


1.0 


0.0138 


0.0183 


40.7 


54.2 



Table 1: Values of Zs, Zp and /„■, calculated with a sharp cut-off, for various 
values of ^ and for two values of M . 

on the shape of the regulator. The ratio Mo/A is therefore the key parameter 
to consider and we must determine the physically meaningful range of values 
for M /A. 

In the vicinity of q = 0, the inverse pion propagator determines the value 
of f n . In the chiral limit, we have: 

Kp 1 (q) « Z P q 2 f v = M ^Z~ P (56) 

and Zp is the residue of the pion propagator at the pole q = 0. For a given 
value of Mp /A, the value of f v depends on the value of M . Tables [I] and 
give the values of Zp and f n for various values of ^ and for two values of 
M , namely M = 300 MeV and M = 400 MeV. For M = 300 MeV, the 
observed value f n = 93 MeV is fitted with « 0.4. For M = 400 MeV, 
higher values ^ ~ 0.6 and ^ ~ 0.7 are required respectively, when a sharp 
cut-off and a gaussian regulator are used. Soliton calculations require M to 
lie between 300 and 400 MeV [[|. Higher values of M and therefore of M /A 
have also been considered in order to push the unphysical qq continuum well 
above the p mass of 770 MeV . In Ref.0 for example, the value M /A = 0.74 
is used. In order to cover the full range of physically meaningful parameters, 
we perform our calculations from M /A = 0.2, which means a relatively high 
value of the cut-off, up to Md/A = 0.8, which means a low value of the 
cut-off. 

One crucial point here is that the relevant range of parameters (typically 
0.4 < M /A < 0.8) corresponds to uncomfortably low values of the cut-off. 
Most field theoretic methods applied to statistical mechanics and particle 
physics have been developed to systems in which M << A. Considerable 
errors can be (and have been) made by applying methods and concepts, bor- 
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Mo 
A 


Z s 


Z P 


U (MeV) 
forM = 300 MeV 


U (MeV) 
forM = 400 MeV 


0.2 


0.0766 


0.1621 


121 


161 


0.4 


0.0305 


0.0950 


92.5 


123 


0.6 


0.0155 


0.0659 


77.0 


103 


0.8 


0.0089 


0.0498 


67.0 


89.3 


1.0 


0.0056 


0.0395 


59.7 


79.6 



Table 2: Values of Zs, Zp and f w , calculated with a gaussian regulator, for 
various values of ^ and for two values of Mq. 

rowed from the study of systems in which M << A, and applied to systems 
in which the cut-off A is of the same order of magnitude as Mq. Our calcu- 
lations focus on several problems which one encounters when the cut-off is 
not much larger than the calculated observables. This is the regime appli- 
cable to low-energy hadronic physics and it is not an artefact of the Nambu 
Jona-Lasinio model. For example, low energy effective theories derived from 
an instanton liquid |7|, f| yield values of M /A of the order of 0.4. 

As an example, consider the inverse cx-meson propagator at q = in the 
chiral limit. It is given by fl43|): 

Kg 1 (q = 0) = 8N c N f M 2 f M = AM 2 Z S (57) 

where /m is the function f^f (<?), defined in ( |B.5|) , and taken at q = 0, where 
it is independent of n and p. This definition of Zs is quite arbitrary except 
for the fact that, for large values of the cut-off, that is, when — > 0, we 
have Zs = Zp and the Nambu Jona-Lasinio model reduces to a linear sigma 
model. The values of Zs are also listed in the tables [1| and 0. We see that, 
with a sharp cut-off, Zs and Zp differ by about 30% in the relevant parameter 
range ^ = 0.4 — 0.8. With a gaussian cut-off the equality Zs = Zp is not 
even approximately obtained. The fact that Zs ^ Zp contradicts most, if 
not all previously reported calculations of the meson propagators, when they 
are derived from an unregularized action with loop integrals subsequently 
regularized. If we had proceeded this way, the loop integrals f$ (q) would 
have been independent of n and p, and the function g M (q) would have become 
independent of q. Instead of the expression ([43]), the meson propagators would 
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have been equal to the usually quoted expressions (in the chiral limit): 



A7 1 (q) = 4N c N f ~ (q 2 + AM 2 ) f M (q) Kp 1 (q) = 4N c N f ^q 2 f M (q) 

(58) 

and Zs would equal Zp. The tables [l] and [2| show that considerable errors 
can be introduced if the regularization is not specified in the action from the 
outset and adhered to. Of course, these errors would be small if the cut-off 
were large, that is, if M /A were very small. However, in the applications of 
the Nambu Jona-Lasinio model to low energy hadronic physics, it is not. 

The recently claimed instability of the Nambu Jona-Lasinio model, her- 
alded by Kleinert et al.jnj, is based on a reduction of the Nambu Jona- 
Lasinio model to a non-linear a-model. Although the arguments presented 
above raise doubts as to the validity of this reduction, and such doubts have 
also been voiced elsewhere ||15|| , we shall show in section [10] that instabilities 
do indeed arise but that they are artefacts of the sharp cut-off regularization 
when the cut-off is close to M . 



8 The quark condensate of the unconstrained 
system. 



In this section we consider the quark condensate (|53| ) of the unconstrained 
system (j = 0). Various contributions to the quark condensate are listed in 
table for various values of Mo/A. They are expressed in units of A 3 . We see 
that, throughout the range of relevant parameters, the quadratic fluctuations 
of the fields do not alter significantly the quark condensate. They show no 
sign of restoring chiral symmetry. These results agree with those found in 



Ref.||12||. A finer analysis would show that the negative contributions of the 
pion field are due to the exchange term. Although the field fluctuations do not 
alter significantly the ground state expectation value of the quark condensate, 
we shall see in the next section that they do cause an appreciable quantum 
fluctuation of the condensate. 
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Mi 

A 


\ 1 class 


\ /a loop 


\ 1 IT loop 


\ 1 total 


0.4 

sharp 


-0.0392 


0.0018 


-0.0030 


-0.0404 


0.6 

oho VT~\ 

blldl p 


-0.0446 


0.0004 


0.0003 


-0.0439 


0.8 

sharp 


-0.0451 


-0.0001 


0.0020 


-0.0432 


0.4 

gauss 


-0.0218 


0.0024 


-0.0026 


-0.0220 


0.6 

gauss 


-0.0275 


0.0016 


-0.0011 


-0.0270 


0.8 

gauss 


-0.0315 


0.0012 


-0.0001 


-0.0304 



Table 3: The classical quark condensate and the contributions of the 
quadratic fluctuations of the a and ir fields. The last column sums all the 
contributions. The condensates are expressed in units of A 3 . 

9 The quantum fluctuations of the quark con- 
densate in the unconstrained system. 

The expressions and ([49]), which give the quark condensate u/>V) an d 

the squared condensate / (ipT a ip^J \ , allow us to calculate the quantum fluc- 
tuation of the quark condensate A^,: 

= ^(^) 2 )-(^) 2 (59) 

The values are listed in table |] for various values of Mo/ A. The fluctuations 
are due to the exchange and ring diagrams and they vanish in the classical 
approximation. We see that, relative to the quark condensate, they are quite 
large: about 50% when a sharp cut-off is used and between 70% and 80% 
when a gaussian regulator is used. We also see that the quark condensates are 
more sensitive to the shape of the regulator than f n . This is because would 
diverge logarithmically with a large cut-off, whereas the quark condensates 
would have a quadratic dependence on the cut-off. This is also the reason 
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Mo 

A 

A 


\ / 


\ \ T ' 1 






0.2 
sharp 


-0.0300 


11.69 x 10~ 4 


0.0164 


0.55 


0.4 

sharp 


-0.0405 


20.05 x 10~ 4 


0.0192 


0.47 


0.6 

sharp 


-0.0438 


23.32 x 10~ 4 


0.0200 


0.46 


0.8 

oho VT~\ 

snai p 


-0.0432 


22.56 x 10~ 4 


0.0196 


0.45 


0.2 

gaussian 


-0.0145 


4.38 x 10- 4 


0.0151 


1.03 


0.4 

gaussian 


-0.0220 


7.75 x 10" 4 


0.017 


0.78 


0.6 

gaussian 


-0.0270 


10.67 x 10~ 4 


0.018 


0.68 


0.8 

gaussian 


-0.0304 


12.97 x 10~ 4 


0.019 


0.63 



Table 4: The quark condensate (ipip), in units of A 3 , the squared condensate 



(^)T a i/j\ V in units of A 6 , the fluctuation of the condensate A, in units of 



A and its relative value A^,/ (ijjijj 



why, once f n is fixed, larger condensates are obtained with a sharp cut-off 
than with a gaussian regulator. 



10 The effective potential in the chiral limit. 

Figure |T] shows the classical effective potential ( |56] ) in the chiral limit, as a 
function of -r, calculated with a gaussian regulator, for various values of 
The minimum occurs at M = Mo and in all this work, we define the zero 
of energy to be equal to the minimum of the classical action at the point 
M = Mq. These curves map out the energy surface of the system while it is 

being deformed by the constraint |j (^pT a tjj} . As stressed in section |5.1| , it 
makes no difference whether we plot the effective potential as a function of 
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j, or of M. Plotted as a function of M, the curves are easier to understand. 
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Figure 1: The classical effective potential calculated with a gaussian regulator 
and plotted as a function of M, for various indicated values of M /A. The 
effective potential is expressed in units of A 4 . 



Figure || shows shows the classical effective potential (|36|) calculated with 
a sharp cut-off. The code can be checked against analytic expressions in this 
case. We see that for increasing values of that is, for decreasing values 
of the cut-off, the minimum of the effective potential at M — M becomes 
increasingly shallower and that it disappears altogether at the critical value 
^ > 0. 742 which can be evaluated analytically. When a gaussian regulator 
is used, the onset of the instability occurs at the much higher value ^ > 2.93 
which was evaluated numerically. 

The system appears to display an instability with respect to perturbations 
caused by the constraint |j (^pT a il?J . Furthermore, the energy of the system 
does not seem bounded from below. It goes without saying that the classical 
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Figure 2: The classical effective potential, calculated with a sharp 4- 
momentum cut-off, and plotted as a function of M for various indicated 
values of Mo/A. The effective potential is expressed in units of A 4 . 

Let us take a closer look at this apparent instability. It is not an artefact 
of the classical approximation. Figure [3] shows the various contributions to 
the effective potential when the l/N c field fluctuations are included, using a 
sharp cut-off with Mo/A = 0.8. We see that the field fluctuations lower the 
energy but that they do not significantly change the shape of the effective 
potential, so that the instability remains. We also see that the exchange 
(Fock) term dominates the 1/N C corrections because it is more sensitive than 
the ring diagrams to the high momenta running in the meson loop. These 
conclusions remain valid for smaller values of Mo/A. 

By way of comparison, figure ^ shows the contributions of the 1/N C cor- 
rections when a gaussian regulator is used. Although the effect of the meson 
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Figure 3: Various contributions to the effective potential calculated with a 
sharp cut-off with a value M /A = 0.8. The energies are measured relative 
to the classical action at M = M . 

fluctuations is somewhat larger, the stability of the system is not modified. 
We see however that the ring diagrams make a somewhat larger contribution 
than the exchange (Fock) term, because the gaussian regulator reduces the 
effect of the high momenta running in the meson loop. The same conclusions 
can be reached for different values of As mentioned above, the instabil- 
ity also occurs when a gaussian cut-off is used, but at much higher values of 
Mo/ A > 2.93. For such high values, the cut-off is too small to be physically 
meaningful. With a gaussian regulator and in the relevant range of parame- 
ters 0.4 < M /A < 0.8, one needs to probe the system with values as high as 
M/A > 4 before it becomes apparent that the energy is not bounded from 
below. 

A clue concerning the nature of the instability can be obtained by con- 
sidering the effective potential obtained with a sharp 3-momentum regular- 
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Figure 4: Various contributions to the effective potential calculated with a 
gaussian regulator and with M /A = 0.8. The energies are measured relative 
to the classical action at M = Mq. 



isation. In this regularization, the trace of the quark loop is calculated by 
integrating the energy variable from — oo to +oo, and by limiting the 3- 
momentum by the condition k < A. This regularization is tantamount 
to a limitation of the quantum mechanical Hilbert space available to the 
quarks. Figure [5] shows that the effective potential, calculated with a sharp 
3-momentum cut-off, behaves as expected and that it does not display the 
instability. 



11 Unphysical poles of the quark propagator. 

The fact that the effective potential, calculated with a sharp 4-momentum 
cut-off, displays an instability that does not occur when a sharp 3-momentum 
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0.010 




Figure 5: The classical effective potential calculated with a sharp 3- 
momentum cut-off, for various indicated values of Mo/A. No instability ap- 
pears. 



cut-off is used, can be understood as an effect of the unphysical poles of the 
quark propagator which are introduced by the 4-momentum regulator. For 
constant fields, the quark propagator can be written in the form: 

1 l —H7 = ~ k ^ + rlM (*„ = U k)) (60) 

k^ + r 2 M u 2 + k 2 + riM 2 ^ \ 1 J) K J 



When a 3-momentum regularisation is used, we have r\ — 1 when 



< A 



and r\ = otherwise. In the complex u;-plane, the quark propagator has only 



on-shell poles at lu = ±iy k 2 + M 2 . 

However, when a 4-momentum regulator is used, the quark propagator 
acquires extra poles, which also occur when proper-time regularization is used 



[ 16| . Such poles are unphysical in the sense that, taken seriously, they lead to 
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—0.01 1 T i 0-04 7 











Table 5: Position and residues of the poles of the quark propagator with a 
gaussian regulator. Values are given for M/A = 0.2 and 0.8 



instabilities of the vacuum. Equivalently one can say that the system behaves 
as if it was governed by a non-hermitian hamiltonian. This is sometimes 
also expressed by saying that the theory becomes acausal. We assign the 
cause of the instability discussed above to the existence of such unphysical 
poles. When the 4- momentum cut-off is high (and Md/A correspondingly 
low) the effect of these unphysical poles is not felt. But this is not the case 
in low energy hadronic physics where the parameter Mo/A is in the range 
0.4 < Mo/A < 0.8. Furthermore, the position (and therefore the effect) of 
the unphysical poles depends very much on the shape of the regulator. 

A qualitative understanding of the difference between a sharp 4-momentum 
cut-off and a soft gaussian regulator can be understood by comparing the lo- 
cation of the poles of the quark propagator in the complex k 2 plane. When 
a gaussian regulator is used, the poles of the quark propagator occur when: 

k 2 + M 2 e~^=0 (61) 

Poles on the real axis of the complex k 2 plane occur only if M 2 < 4j-. Oth- 
erwise and in addition, poles occur in the complex plane. Figure ^ shows 
the poles of the quark propagator in the two cases M/A = 0.2 (large cut-off) 
and M/A = 0.8 (small cut-off). The poles all lie to the left of the imaginary 
axis and they move closer to it when the cut-off gets small. Table [| gives the 
residues of the poles. 

Consider next regularisation using a sharp 4-momentum cut-off. The cor- 
responding regulator does not have an analytic form, but we have checked 
that very similar results are obtained with a Wood-Saxon shaped regulator 
r t = l^-ffi which becomes equivalent to a sharp cut-off when c — > 0. The 

l+e c 
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Figure 6: The poles of the quark propagator with a gaussian regulator in the 
complex k 2 plane. The circles are for M/A = 0.2 and the diamonds are for 
Mj A = 0.8. 



poles of the quark propagator are then the solution of the the equation: 

M 2 



k z + 



1 + e" 



fe2_ A 2 







(62) 



Setting k 2 = x + iy, this equation decomposes into the two equations: 



x + M z 



y-M 2 



1 + e c cos - 



1 + e c cos * ) + I e c sin * 



e c sin 



in ^ 



x-A2 



1 + e c cos - ) + e <= sin 



-A 2 



(63) 



(64) 
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y 


residue 


M = 0.8 c=0.05 


-0.64000 
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M = 0.8 c=0.05 


1.02152 


±0.46354 


-0.02 7646 ± 0.007 3992i 


M = 0.8 c = 0.05 


1.01790 


±0.77472 


-0.02 4736 ±0.01 1194z 


M = 0.8 c = 0.05 


1.01431 


±1.08764 


-0.02 1203 ± 0.01 3627i 


M = 0.8 c = 0.05 


1.01128 


±1.401618 


-0.01 7739 ± 0.0 U8i 



Table 6: Position and residues of the first few poles the quark propagator 
when a Wood-Saxon shaped regulator with c = 0.05 A 2 and M/A = 0.8. 



One pole occurs on the negative real axis (y — x « — M 2 ) in the vicinity 
of — M 2 . It is simple to see that the complex poles all occur in the vicinity 
of the line x — A. The imaginary part y is then the solution of the equation: 

— = tan — (65) 

2c 4c 2c y ' 

The spacing between the poles is thus close to 27rc so that they get denser 
in number as c —>■ 0. In that limit, the continuum of poles forms a cut which 
expresses the discontinuity of the regulator when a sharp cut-off is used. A 
more exact numerical calculation of the position and residues of the poles is 
given in the table [I] for the case where c = 0.05 A 2 and M/A = 0.8. 

The unphysical poles produced by a soft gaussian regulator lie to the left 
of the imaginary axis of the complex k 2 plane. They are therefore mostly 
felt at low values of k 2 where phase space factors reduce their effect. The 
unphysical poles produced by a sharp cut-off lie close to the boundary k 2 = A 2 
of high values of k 2 from which diverging quantities derive most of their 
contribution. This explains qualitatively the difference between the effect of 
the two regularizations on the effective potential. It would be worth analyzing 
whether the instabilities, recently heralded by Kleinert et al. [?], are not also 
artefacts of the use of a sharp 4- momentum in conjunction with a low cut-off. 



12 Conclusions. 

Although the meson loop contributions do not modify appreciably the value 
of the quark condensate, they do cause large quantum fluctuations of the 
quark condensate. The Lorentz invariant regularization of the quark propa- 
gator, used in conjunction with the relatively low cut-off values required in 
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low energy hadronic physics, makes the physical vacuum unstable against 
distortions caused by the squared quark condensate. The instability depends 
strongly on the shape of the regulator. It is only weakly felt when a soft 
gaussian regulator is used but its effects are greatly enhanced in calcula- 
tions which use a sharp 4-momentum cut-off. Large errors can be made if 
loop integrals are regularized after being derived from an unregularized ac- 
tion, instead of including the regulator in the model action from the outset. 
The ground state instability can be traced to unphysical poles of the quark 
propagator which are introduced by the regulator. 

A Expressions for the classical action and the 
classical effective action. 

In the chiral limit, the value of the classical action at the stationaty point is: 

I hm =o (M 2 ) - / (Ml) (A.l) 

= 2 c f (2vr) 4 J 4 v k 2 + r\M 2 + k 2 + r\M 2 k 2 + r\M 2 ) 

(A.2) 

We measure all energies relative to the minimum / (Mg ) of the classical action 
in the unconstrained system where j = 0. In units of f2A 4 , the classical action 
depends only on the two variables ^ and 



B Expressions for the meson propagators. 

B.l The meson propagators and the polarization 
function Il a at a point ip a = (M, 0, 0, 0). 

From the second order expansion of —Tr In (—id fl '~f fl + r(p a T a r) we obtain the 
following expression for the polarization function: 
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The calculation is standard. Taking traces over the Dirac and flavor indices, 
and keeping track of the regulator, we obtain: 



n( 2 ) = 4ivyv7± £ 5S {q) 6S (- q ) (lff% ( q ) + M 2 (f™ ( q ) + ( q )) - 9M ( q ) 
+ 4A ^/^ E ^ (?) ^ (-«) & 2 /£ (?) + m2 (/m (?) - /m (?)) - <?M (?) 



(B.2) 

where we wrote the fields f a in terms of scalar and pseudoscalar fields S and 

T a ip a = S + i^nPi (B.3) 
The Fourier transforms are defined to be: 

S{q) = J d 4 x e iq ^S (x) Pi (q) = J d 4 x e iq ^P { (x) (B.4) 

The function f^f (q) is: 

Q(q) = j k 3 dk dasin 2 a 7 ^ r (B.5) 

JMKq> (27r) 4 io Jo (kl + rtM^)(kl + ri 2 M^) 1 ' 

with ki — k — | and k 2 — k + |. The function g M (q) is: 

^.yp roc rir 

9M (q) = J^J Q k'dk J o da sin a -, - (B.6) 

and we denote by qm the function gM (q — 0). 

It follows that the polarization function is diagonal in momentum space: 

(qa\U\q'b) = 5 ab 5 qq ,U a (q) (B.7) 

where: 

n a=0 (q) = U S (q) = AN c N f {^-q 2 f 2 (q) + M 2 (/* (q) + f% (q)) - g M (q) 

n a=1;2 , 3 (q) = U P (q) = 4N c Nf ^-q 2 ^ (q) + M 2 (/* (q) - (qj) - g M (q) 

(B.8) 
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The inverse propagator matrix K 1 is obtained by adding — (xa (V — j) 1 ybj. 
We use the gap equation ( p?D to write: 



~ (V - jy 1 = W c N fJ ^L-g M (B.9) 



so that: 



Ks 1 (q) = 4N c N f (\ff$ (q) + M 2 {f* (q) + f$ (q)) - g M (q) + ^^Qu (0) 



M -m- 

K? (q) = AN c Nf Qg 2 / 22 (?) + M 2 (g) - f$ (q)) - g M (q) + j^9m (0) 

(B.10) 
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